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We develop a unified, information theoretic interpretation of the number-phase complementarity 
that is applicable both to finite-dimensional (atomic) and infinite-dimensional (oscillator) systems, 
with number treated as a discrete Hermitian observable and phase as a continuous positive oper- 
ator valued measure (POVM). The relevant uncertainty principle is obtained as a lower bound on 
. . . entropy excess, X, the difference between the entropy of one variable, typically the number, and 

. the knowledge of its complementary variable, typically the phase, where knowledge of a variable is 

■ defined as its relative entropy with respect to the uniform distribution. In the case of finite dimen- 
' sional systems, a weighting of phase knowledge by a factor fi (> 1) is necessary in order to make 

the bound tight, essentially on account of the POVM nature of phase as defined here. Numerical 
and analytical evidence suggests that fi tends to 1 as system dimension becomes infinite. We study 

' the effect of non-dissipative and dissipative noise on these complementary variables for oscillator as 

^ well as atomic systems. 

■ PACS numbers: 03.65.Ta,03.65.Yz,03.67.-a 

I. INTRODUCTION 

> ■ 

^ i Two observables A and i? of a d-level system are called complementary if knowledge of the measured value of A 
^ ' implies maximal uncertainty of the measured value of B, and vice versa Complementarity is an aspect of the 
Q^l Heisenberg uncertainty principle, which says that for any state -0, the probability distributions obtained by measuring 
A and B cannot both be arbitrarily peaked if A and B are sufficiently non-commuting. Expressed in terms of 
I measurement entropy the Heisenberg uncertainty principle takes the form: 

H{A) + H{B)> log d. (1) 

, where H{A) and H{B) are the Shannon entropy of the measurement outcomes of a d-level quantum system 

ff^ ' Eq. ([1]) has several advantages over the traditional uncertainty multiplicative form [^-0] . 

jy-^ , More generally, given two observables A = a|a)(a| and B = h\h){h\, let the entropy generated by measuring 

■ A or B on a state 1-0) be given by, respectively, H{A) and H{B). The information theoretic representation of the 

. Heisenberg uncertainty principle states that H{A)+H{B) > 2 log (^jjXb)) ' '^here f{A, B) = maxo^h |(a|6)|, and ff (•) 

• • ■ is the Shannon binary ent ropy . A pair of observables, A and B, for which f{A, B) — d~^/^ are said to form mutually 
. ^ ' unbiased bases (MUB) [13] • Conventionally, two Hermitian observables arc called complementary only if they are 
, mutually unbiased. 

An application of this idea to obtain an entropic uncertainty relation for oscillator systems in the Pegg-Barnctt 
5^ , scheme Q has been made in Ref. [ll|, and for higher entropic uncertainty relations in Ref. ^3- An algebraic 
treatment of the uncertainty relations, in terms of complementary subalgebras, is studied in Ref. [IJ 



An extension of Eq. ([T]) to the case where A or B is not discrete is considered in Ref. [ij] , where the problem that the 
Shannon entropy of a continuous random variable may be negative is circumvented by instead using relative entropy 
(also called Kullback-Leibler divergence, which is always positive) [isl ITg} with respect to a uniform distribution. 
This quantity is a measure of knowledge [l^l ■ An example of where this finds application would be when one of the 
observables, say A, is bounded, and its conjugate B is described not as a Hermitian operator but as a continuous- 
valued POVM. A particular case of this kind, considered in detail in Ref. is the number and phase of an 
atomic system. This generalization of the entropic uncertainty principle to cover discrete-continuous systems still 
suffers from the restriction that the system must be finite dimensional, since in the case of an infinite-dimensional 
system, such as an oscillator, entropic knowledge of the number distribution can diverge, making it unsuitable for 



'Electronic address: 'subhashish@cm i. acini 
t Electronic address: ,srik@ppisr.res.in] 



2 



infinite-dimensional systems. Therefore to set up an entropic version of the uncertainty principle, that unifies and is 
applicable to all systems, including infinite dimensional and/or continuous- variable systems, it may be advantageous 
to use a combination of entropy and knowledge, in particular, the difference between entropy of the discrete, infinite 
observable and between phase knowledge. This is discussed in detail below. 

The theory of open quantum systems addresses the problems of damping and dephasing in quantum systems by its 
assertion that all real systems of interest are in fact 'open' systems, each surrounded by its environment. One of the 
first testing grounds for open system ideas was in quantum optics [17j . Depending upon the system- reservoir (S — R) 
interaction, open systems can be broadly classified into two categories, viz., quantum non-demolition (QND), where 
[Hs,Hsr] = resulting in pure decoherence, or dissipative, where [Hs,Hsr] ^ resulting in decoherence along with 
dissipation [isj . 

The plan of the paper is as follows. In Section II, we briefly introduce, in anticipation of the discussion to follow, the 
concept of quantum phase distributions for oscillator as well as two-level atomic systems. In Section III, we develop 
an information theoretic representation of complementarity. A central feature here is the study of number-phase 
complementarity using the principle concept of entropy excess, the difference between number entropy and phase 
knowledge, mentioned above. The use of the entropy excess enables a unifled, information theoretic interpretation 
of the number-phase complementarity, with dimension-independent lower bound, that is applicable both to flnite- 
dimensional (atomic) and infinite-dimensional (oscillator) systems, as well as discrete (number) and continuous (phase) 
variables. 

We apply this entropic uncertainty principle to various physical systems: oscillator systems (both harmonic as well 
as anharmonic), in Section IV, and atomic systems, in Section V, for a host of physically relevant initial conditions. In 
addition, the effect of purely dephasing as well as dissipative influences on the system's evolution, due to interaction 
with its environment, and hence the entropy excess is studied for each case considered in Sections IV and V. In Section 
VI, we make our conclusions. 



II. QUANTUM PHASE DISTRIBUTIONS 

The quantum description of phases [l^ has a long history [sl. [20l - [24j : see also Refs. [2^, [2^. In a recent approach, 
which we adopt, the concept of phase distribution for the quantum phase has been introduced (25l.[27j. Here we briefly 
recapitulate, for convenience, some useful formulas of quantum phase distributions for oscillator systems (28l. f29l|. For 
the case of atomic systems, the basic formulas were presented in [T^ . 

Following Agarwal et al. j27| we define a phase distribution 7^(0) for a given density operator p, which in our case 
would be the reduced density matrix, as 



r £ Pn^,ne'("-™>^ (2) 



27r 

m,n— 

where the states |^) are the eigenstates of the Susskind-Glogower [2lj phase operator corresponding to eigenvalues of 
unit magnitude and are defined in terms of the number states \n) as 

oo 

|0)==^e™^|n), (3) 

n=0 

The sum in Eq. ([2]) is assumed to converge. The phase distribution is positive definite and normalized to unity with 

!e\()){o\de = i. 

The complementary number distribution is 

p(to) = (m|p|TO), (4) 



where |m) is the number (Fock) state. Analogous results exist for atomic states, with the Susskind-Glagower states 
replaced by atomic coherent states [sil [s!] , and number states by Wigner-Dicke states [11] . 
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III. INFORMATION THEORETIC REPRESENTATION OF COMPLEMENTARITY 

Defining entropic knowledge R[f] of random variable / as its relative entropy with respect to the uniform distribution 



i le 



R[f] ^ S = ^/(j)log(d/0-)), (5) 

we can recast Heisenberg uncertainty principle in terms of entropy H and knowledge i?, as shown by this easy theorem 

Theorem 1 Given two Hermitian observables A and B that form a pair of MUB, the uncertainty relation can be 
expressed as 

X{A,B) = H{A)-R{B)>0. (6) 

Proof. Let the distribution obtained by measuring A and i? on a given state be, respectively, {pj} and {qt}- Denoting 
H{A) = J2j Pj log2Pj' the l.h.s of Eq. © is given by 

H{A)-s(^B\\^^ = HiA)-J2qklog{dqk) 

= H{A) + H{B)-\ogd (7) 
> 21ogY ^, )-\ogd. (8) 



\fiA,B) 

where Eq. ^ follows from Ref. For a pair of MUB d, f{A, B) ~ d^^/^, from which the theorem follows. ■ 

From Eq. ([7]) it follows that X{A,B) = X{B,A). Therefore, phyically Eq. ([6]) expresses that ignorance of one of 
two MUB variables is at least as large as the knowledge of the other. It is not difficult to see that X{A,B) attains 
its largest value of log(o?) when A and B are MUBs, and its minimum value of — log((i) when A and B are identical. 
This gives a way to quantify the 'degree of complementarity'. Define Xnun{A, B) are the smallest value of X{A,B) 
over all possible states for a given pair of Hermitian observables A and B. Then, two observables A and B are 
maximally complementary (i.e., MUB) if Xy^i^{A, B) ~ 0, and they arc minimally complementary (i.e., identical) if 
X,,i„(AS) = -log(d). 

A point worth noting about Eq. ([6]) is that it contains no explicit mention of dimension d. What is remarkable 
is that we find this situation persists even when one of A or B is not discrete, but a continuous-valued POVM 
(for discrete- valued POVMs. cf. Rcf. [1^), and furthermore, the system is no longer finite dimensional but instead 
infinite dimensional. The only additional requirement is that the continuous-valued variable should be set as B (the 
knowledge- rather than the ignorance- variable) , since H{B) can potentially be negative for such variables. This 
makes X{A,B) > as a very succinct and general statement of the uncertainty principle. By contrast, because 
there is no prior guarantee that measurement entropy H{-) will be non- negative for a continuous- valued observable, 
it is not obvious that the version of the Heisenberg uncertainty principle given by ([1} is generally applicable, and 
furthermore, because there is no prior guarantee that measurement entropic knowledge R{-) will be well-defined for 
infinite-dimensional variables, the version R[A) + R{B) < \og{d) of Ref. is also not obviously generally applicable. 

One catch is that on account of the POVM-nature of B, R{B) may have a maximum value less than log(rf) in the 
finite dimensional case. It will be to generalize the concept of 'maximal complementarity' or 'MUBness' to apply 
those terms to A and B, when one of them is a POVM, if the maximal knowledge of the measured value of A implies 
minimal knowledge of the measured value of B, and vice versa, but with maximum knowledge no longer being required 
to as high as log d bits. 

For the phase variable given by the POVM cj) and probability distribution 'P{(j>), entropic knowledge is given by the 
functional [Hiil: 

R[Vm= dct>Vmog[2T:V{^)l (9) 



where the log(-) refers to the binary base. 
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It is at first not obvious that Eq. ([6]) holds for infinite dimensional systems. Based on a result due to Ref. [S^ for 
an oscillator system, which in turn uses the concept of the {p, g)-norm of the Fourier transformation found by Beckner 
[soj for all values of p, for an oscillator system, we can show that it is indeed the case. In particular, 

- r d<t>P{4>) \0g{Pm -f^Pn^ \0g{p,n) > l0g(2^) (10) 

Setting the 'number variable' m in Eq. (|10p as A, and the phase variable (j) as B, and noting that the first term in 
the l.h.s of Eq. dTU]), using Eq. ©, is just log(27r) - i?[P(0)], we obtain 

X[m, = H[m] - R[(t)] > 0, (11) 

which is Eq. © applied to an infinite-dimensional system that includes a non-Hermitian POVM (phase (j)). Eq. 
(|lip expresses the fact ignorance of variable m is at least as great as knowledge of its complementary partner, (j). 
Comparing Eqs. ^ and pT|). we find that the statement X > as a description of the Heisenberg uncertainty 
relation holds good both for finite and infinite dimensional systems. The version X > of the Heisenberg uncertainty 
principle may be called the principle of entropy excess. An information theoretic interpretation of the above relation 
has been studied, in the context of phase resolution in harmonic oscillator systems, in [40j . Also, the number-phase 
complementarity, for a harmonic oscillator system, using information exclusion relations has been studied in [4l| . 



IV. OSCILLATOR SYSTEM 



Here we consider the application of the principle of entropy excess (|lip to oscillator systems, both harmonic as well 
as anharmonic, starting from a number of physically relevant and interesting initial conditions and interacting with 
their environment via a purely dephasing (QND) as well as dissipative interaction. The strategy would be to compute 
the phase and number distributions for each case, use them to obtain phase knowledge ([9|), number entropy and use 
them in Eq. (jlip to study the entropy excess and thus the number-phase complementarity in oscillator systems. 

A. QND system-bath interaction 

Consider the following Hamiltonian describing the interaction of a system with its environment, modelled as a 
reservoir of harmonic oscillators, via a QND type of coupling : 

H = Hs + -ffif + HsR 

Here Hs, H]^ and Hsr stand for the Hamiltonians of the system, reservoir and system- reservoir interaction, respec- 
tively. Hs is a generic system Hamiltonian which we will specify in the subsequent sections to model different physical 
situations. bk denote the creation and annihilation operators for the reservoir oscillator of frequency Wfe, gk stands 
for the coupling constant (assumed real) for the interaction of the oscillator field with the system. The last term on 
the right-hand side of Eq. (1) is a renormalization inducing 'counter term'. Since [Hs,Hsb] = 0, the Hamiltonian 
(1) is of QND type. The system plus reservoir composite is closed obeying a unitary evolution given by 

Pit) = e-T."'p{Q)ei"\ (13) 

where 

p{0)^p'{Q)pr{0), (14) 

i.e., we assume separable initial conditions. The reservoir is assumed to be initially in a squeezed thermal state, i.e., 
it is a squeezed thermal bath, with an initial density matrix pr(0) given by 

PRiO) ^ S{r,<i>)ptJHr,'f), (15) 

where 

Pth = n [1 ~ e-'"*""] e-'^^'^^'^i'^" (16) 
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is the density matrix of the thermal bath, and 



in 



(17) 



is the squeezing operator with r^, $fc being the squeezing parameters j42l |. We are interested in the reduced dynamics 
of the 'open' system of interest S, which is obtained by tracing over the bath degrees of freedom. Using Eqs. (fT2l). (fT4|) 
in Eg. (jl3p and tracing over the bath variables, we obtain the reduced density matrix for 5, in the system eigenbasis, 
as [131 



'lit) PnrM- 



(18) 



In the above equation. En is the eigenvalue of the system in the system eigenbasis while r]{t) and j(t) quantify the 
effect of the bath on the system and are given in Appendix El for convenience. 



1. System of a harmonic oscillator 



We consider the system 5 of a harmonic oscillator with the Hamiltonian 

Hs = huj [a^a + . (19) 

The number states serve as an appropriate basis for the system Hamiltonian and the system energy eigenvalue in this 
basis is 

En = nw (^i + . (20) 

The harmonic oscillator system is assumed to start from the following physically interesting initial states: 
(A). System initially in a coherent state: 
The initial density matrix of the system is 

p^(0) = |«)(a|, (21) 

where 

a = \a\e'^'> (22) 
is a coherent state Making use of Eqs. (fT8)) . ((2T|) in Eq. ([2|), the phase distribution is obtained as [28j 

In+m 



-e ^-i(m-n)(e-eQ) ^-iui(m-n)t 



27r ^ \Jn\m\ 

m,n— 

^^i(huif (m-n)(n+m+l)i)(t) ^-{huj)'^ (n-mf -i(t) ^ (|23) 

The corresponding complementary number distribution is obtained, using Eq. as 

lap™ I i2 

p(to) = ^— e-l"l . (24) 
m! 

Using V{9) in Eq. ^ to get the phase knowledge, p{ra) ((M)) to get the number entropy and using these in Eq. 
(|TT1) we get the entropy excess. These are plotted in Figures [T] It is clearly seen, by a comparison of Figure [TJb) with 
(a) (representing unitary evolution), that including the environmental effects due to finite temperature and squeezing 
causes the entropy excess to increase by randomizing phase and thus causing R[d] to fall, whereas H\rn\ remains 
invariant because QND interactions characteristically leave the number distribution pirn) Q invariant [2^ . This can 
be seen from Eq. (j24p . where the only parameter entering the distribution p(rn) is the initial state parameter a. The 
figures clearly show that the principle of entropy excess, Eq. (|lip. is satisfied for both unitary evolution as well as in 
the case of interaction with the bath. 

(B). System initially in a squeezed coherent state: 
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FIG. 1: Number entropy _ff[?n] (large-dashed line), phase knowledge R[9] (small-dashed line) and entropy excess X[m, ^] (Eq. 

bold line) plotted as a function of the parameter a (|22p for the harmonic oscillator system initially in a coherent state. 
Figure (a) represents the case of the pure state case. We note that as number increases, with increase in a, so does H[m] (since 
the variance of a Poisson distribution equals its mean), whereas phase becomes increasingly certain, leading to increase in 
R[(f>]. Figure (b) represents the case of the system subjected to QND interaction with the parameters lj — 1.0, lJc — 100, 70 
(|A1|) = 0.0025, |ap = 5, = (|22p and with bath squeezing parameters (|A5|I r = 2.0 and a = for a temperature T (in units 
where h = ks = 1) ~1 and an evolution time t = 0.5. 



The initial density matrix of the system is 

p^(0) = |^,a)(a,^|, (25) 

where the squeezed coherent state is defined as [1^ 

\^,a)^S{ODia)\0). (26) 

Here 5* denotes the standard squeezing operator with ^ = rie''^ and D denotes the standard displacement operator 
[iaj . Making use of Eqs. dUl), (gS]) in Eq. (g]), the phase distribution is obtained as ^ 



(tanh(ri)) * 2 ' 



m.„=o 2^-^^/^\ cosh(ri) 



X exp [-|ap(l - tanh(ri) 008(26*0 - 1/'))] 



•\/sinh(2ri) 



Y^sinh(2ri) 



^ ^—iuj{m — n)t^i{h(jj)' (m — n) {n~\-m~\-l)ri{t) ^ — {huj) (n — m)"7(£) 



(27) 

Here i?„[z] is a Hermite polynomial. The corresponding complementary number distribution is obtained, using Eq. 
(H, as 



p{m) 



1 (tanh(ri))' 
2™m! cosh(ri) 



■ exp [-|ap(l - tanh(ri) cos(26'o - -0))] \H„ 



|a|e'(^°-T) 
■\/ sinh(27'i) 



(28) 



Using V{0) (P7)) in Eq. to get the phase knowledge, p{m) ([25)1 to get the number entropy and using these in Eq. 
([TT|) we get the entropy excess which are plotted in Figures [5] From the Figures [2] it can be seen that phase gets 
randomized, resulting in a fall in the phase knowledge R[0], with increase in the system squeezing parameter ri (j26p . 
The number entropy H[m] is not effected by the reservoir, due to the QND nature of the interaction but as can be 
seen from Eq. ([28]), the number distribution p(m) depends upon the initial state parameters a, ri and ip. Thus H[m\ 
as a function of the system squeezing parameter ri first falls and then rises as a result of which the entropy excess at 
first goes down and then rises. The principle of entropy excess, Ec^. is clearly seen to be satisfied. 

An interesting feature here is that in Figure Eljb), even though in comparison with the settings in Figure 
temperature T has increased, the value of R[9] has also increased, contrary to the expectation that temperature 
would cause phase to randomize and thus reduce R[9]. The reason is that the P{9) distribution at T = has a 
bimodal (double-peaked or double-bunched) form, having relatively large variance and thus low R[9]. As temperature 
is increased to T = 1, this bimodal distribution at first collapses into a single-peaked form, the resulting sharp 
reduction in variance, being responsible for the rise in R[9]. With further increase in temperature, the expected 
diffusion of the phase sets in, and R[9] registers a gradual reduction. 
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(a) (b) 

FIG. 2: Number entropy _ff[?n] (large-dashed line), phase knowledge R[9] (small-dashed line) and entropy excess X[m, ^] (Eq. 
((TT}) (bold line) plotted as a function of the system squeezing parameter ri (|26|l for a harmonic oscillator system initially in a 
squeezed coherent state and subjected to a QND interaction, with parameters ui = 1, ujc = 100, 70 = 0.025, lap = 5, = 
(|22|l . ^ (|26p — and with bath squeezing parameters HA5|) r = 1.0 and a = 0. Figure (a) represents an evolution time t — 0.1 
and T = 0, while figure (b) depicts the case for an evolution time t = 0.1 and T = 1. 



2. System of an anharmomc oscillator 



We consider the system S of an anharmonic oscillator with the Hamiltonian 



i/5 = fi^(ata+i)+^(at)V. (29) 



As shown in [28[, the above Hamiltonian can be expressed in terms of the generators of the group SU(1,1) as a 

i or ^ 



result of which the appropriate basis for it would be |m, k) where m = 0, 1, 2, ... and k equal to \ or |. The case of 



k = J corresponds to states with even photon number with the vacuum state coinciding with the vacuum state of the 
harmonic oscillator, while the case of fc = | corresponds to states with odd photon number. Using the properties of 
the SU(1,1) group generators, the action of Hs ((29)) on the basis is found to be 

Hs\m,k) = 2h[uj{m + k) + \m{m + 2k-l)]\m,k) 

= E,njm,k). (30) 

We make use of this to obtain the phase distribution of the anharmonic oscillator system, interacting with a squeezed 
thermal bath via a QND system-bath interaction, and starting from the following physically interesting initial states: 

(A). System initially in a Kerr state: 

The initial density matrix of the system is [3] 

p%0) = \^|JK){M■ (31) 
Here \^jk) is defined in terms of the number states as 

l^if) =^gn|n), (32) 

n 

where 



g„ = ^e^e-^'^"^"-!^ (33) 
Vn! 

In the above equations, |n) represents the usual number state and x = where A is as in Eq. (P^]) . L is the length 
of the medium and v is the speed of light in the Kerr medium in which the interaction has taken place. Making use 



of Eqs. ([T8|. (|3T|) in Eq. ([2]), the phase distribution is obtained as [28 1 



-. 00 

/.TT ' ^ 
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m,n— 
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FIG. 3: Number entropy _ff[?n] (large-dashed line), phase knowledge R[9] (small-dashed line) and entropy excess X[m, ^] (Eq. 

(bold line) plotted as a function of the parameter x (|33|) for an anharmonic oscillator system initially in a Kerr state and 
subjected to a QND interaction. The parameters taken are oj — 1, lJc = 100, 70 = 0.0025, |ap = 5, 6^0 = 0, A = 0.02 and with 
bath squeezing parameters (|A5|I r = 2.0 and a — for an evolution time t = 0.1 and T = 0. 



^ ^4ih^ (m-n)[w+X{m+n- ^)][i^(n+m+^) + X{n^ ~ ^ (m+n))]ri{t) 
^g-4?i='(m-n)2[a;+A(m+n-i)]%(t) 

i_L \^ ,ir7* z2(n-m)^ -22(m-n)[a;+A(m+n+i)]t 

m.n— 

I )-|-A(n2 -l-m^-l-i (m-l-n))] t7(t) 

^g-4?i^(Tn-ri)^[w+A(in-|-n-|-i)]S(t) 

The corresponding complementary number distribution is obtained, using Eq. as 



p(m) = |g2mP + 192 



m-|-l| 



(34) 



(35) 



where q2rm q2m+i can be obtained from Eq. 

Using 7'(6') ((34| in Eq. ([9]) to get the phase knowledge, p{m) ((35|) to get the number entropy and using these in 
Eq. PT|) we get the entropy excess which are plotted in Figure |31 From the figure , it is evident that as the Kerr 
parameter x increases, phase gets randomized leading to the fall in R[9], whereby entropy excess X[to,0] increases, 
since H[m\ remains unchanged. The invariance of H[m\ under change in the parameter x can be easily seen by using 
Eq. (p3|) in Eq. ((35|) . The principle of entropy excess, Eq. (|TT|) . is clearly seen to be satisfied. 

(B). System initially in a squeezed Kerr state: 

The initial density matrix of the system is 

p'{0) = H'sk){^skI (36) 
Here \4^sk) is defined in terms of the number states as 

\iPsk) =^Sn\n), (37) 



where 



and 



S2n 



P 



q2pG2 



S2m+1 — g2p-|-lG'2m-|-12p-|-l(^), 



(38) 



(39) 



with z = rie^"^ , and Grnp{z) = {m\S{z)\p) , where S{z) is the usual squeezing operator, is given by [4f 



p\m\ \ cosh(ri) y 
^ |^ tanh(ri ^^^"+^^ 



""'2' (sinh(ri))2 



(40) 
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FIG. 4: Plot of number entropy H[m] (large-dashed line), phase knowledge R[9] (small-dashed line) and entropy excess X[m,6] 
(bold line) for an anharmonic oscillator initially in squeezed Kerr state, with respect system squeezing ri p9|l . Figure (a) 
represents the pure state case, while Figure (b) represents the system subjected to QND interaction with a squeezed thermal 
bath with temperature T = 1 and evolution time t = 1. The parameters used are to — 1, tOc = 100, 70 = 0.1, ip = 7r/4, \af = 5, 
9o — 0, X = 0.02, A — 0.02, and with bath squeezing parameters r = 0.1, a = 0. 



Similarly, G2rn+i2p+i{z) is given by 

G2m+12p+l 



(-l)P f {2p + iy.{2m + ly. 



p\m\ 



cosh (ri ) 



^ tanh(ri )V'"+^) 



exp {i{m — p)4') 
3 1 



(sinh(ri))2 



(41) 



Here is the Gauss hypergeometric function [4y|. Making use of Eqs. ([T8| . p6p in Eq. ([2]), the phase distribution 
is obtained as [23 



no) 



^ 00 

- y 

271 ^ 



S2mS2n^ 



i2{n-m)e -2i(m-n)[ui+X{m+n~ ^)]t 



m.n—0 
-. 00 

I ^ A2(n-m)e-2i{m-n)\uj+\lm+n+^}]t 

' ^ ^2m+1^2n+l^ ^ ' 

ra.n—0 

|)+A(ri^+m^ + i(m-Kn))]i)(t) 

^g-4ri^(m-n)^[a;+A{m+n+i)]^7(t) 



(42) 



The corresponding complementary number distribution is obtained, using Eq. (j?]), as 

p[m) = \s2m? + \s2m+l?, 



(43) 

where S2m, S2m+i can be obtained from Eqs. psp and (j39p . respectively. 

Using V{9) dH]) in Eq. (0) to get the phase knowledge, p[m) (^5]) to get the number entropy and using these in Eq. 
(jlip we get the entropy excess which are plotted in Figures |4l In Figure |4] (a), depicting unitary evolution, it is seen 
that phase knowledge almost exactly compensates for the growth of ignorance of number, as a functions of ri, whereas 
, in Figure|4] (b), phase knowledge is rapidly lost, depicting clearly the influence of the environment. The principle of 
entropy excess, Eq. (|lip . is clearly seen to be satisfied for unitary evolution as well as when the anharmonic system 
is interacting with its environment. 



B. Dissipative system-bath interaction 



Here the system-reservoir interaction is such that [Hs, Hsr] ^ resulting in decoherence along with dissipation. 



10 



(A). System of harmonic oscillator interacting with a thermal bath resulting in a Lindblad evolution: 
The initial state of the system is a superposition of coherent states which are 180° out of phase with respect to 
each other [i?} . 

1^) =Ai/2(|a)+e'*|-a)), (44) 

where a. = |a|e"^° and 

A = i[l + cos((/))e-2|"l']-i. (45) 

The state j?/') for = would be an even coherent state and for = tt would be an odd coherent state. The reduced 
density matrix can be shown to have the following form psl ] : 

oo 

p(i)- (46) 



n,m— 



where 



Pn,m(t) = + l (iV(t) + lj 

X E(i-]vw+Tj 4r(i + (-ir+™ + (-i)'[(-i)v"^ + (-ire-l) 

X Fl [-m, / + 1; \N(t){N{t) + l)(sinh(7ot/2))2] . (47) 



Here F'^ is the Gauss hypergeometric function [46|, 70 is a parameter which depends upon the system-reservoir 
coupling strength, 



Qn-^e- — , (48) 



and, 

N{t) = Nth{l - e-''«*), Nth = (e^ " l) ' ■ (49) 

The phase distribution is given by 

^ 00 

^(^) = ^ E (50) 



27r 

m,n— 

where pm,n can be obtained from Eq. (|47| . 

The corresponding complementary number distribution is obtained, using Eq. Q, as 

p{m) = pjnjn{t), (51) 

where pm,m is as in Eq. (|47p . 

Using 7^(6') ([50]) in Eq. ([9|) to get the phase knowledge, p{m) ((5T|) to get the number entropy and using these in Eq. 
(jlip we get the entropy excess which arc plotted in Figures [5l From Figure [5] (a) , pertaining to unitary evolution, 
we note that in the even cat (coherent) state ((/> = 0), ignorance of number approximately equals phase knowledge, 
whereas in the odd cat (coherent) state {<j> = tt) the former significantly outweighs the latter. This thus provides 
a complementaristic characterization of the even and odd cat states. The Figure [5] (b) shows that the effect of the 
dissipative environment causes phase to become randomized, leading to an increased entropy excess at all </> (|44p. The 
principle of entropy excess, Eq. (jlip . is clearly seen to be satisfied, for both unitary as well as dissipative evolution. 

(B). System of anharmonic oscillator weakly interacting with a thermal bath: 

The total Hamiltonian depicting a third-order non- linear oscillator coupled to a reservoir of oscillators [i^ , assumed 
to be initially in a thermal state, is 



H 



(52) 
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FIG. 5: Plot of number entropy H[m] (large-dashed line), phase knowledge R[9] (small-dashed line) and entropy excess X[m, 6] 
(bold line) for a harmonic oscillator, initially in a coherent state superposition H44|l . as a function of the state parameter </!> 
(|44|) . Figure (a) pertains to the pure state case. Figure (b) represents the system subjected to a dissipative interaction with 
the environment for an evolution time t = 0.1 and temperature T = 2. The parameters used are cj = 1, 70 = 0.025, \af = 2, 
(t>o (iMl) = 0. 



The reduced density matrix of the anharmonic oscillator, starting from the initial coherent state |^(0)) = ||^(0)|e"'^''), 
can be solved and made use of to obtain the phase distribution 



i{n — rn)9 



m.n— 

- 00 

= - ^ (m!n!)i/V„,™(i)e 



i{n—m)6 



(53) 



Here 



00 ^ 

fraAt) = exp {[-2iK{m -n) + ^]t) (i?,„_„(t))"+"+i ^ - 



1=0 



Nth + 1 



th 



9rn — n (0 



{m + l)\{n + l)\ 



m, -n; I + 1; (smh(7oAt/2)) 



(54) 



In the above equation, fm+i,n+i{^) contains information about the initial state of the system and for the initial 
coherent state |^(0)) = \\i{Q)\e"^°) is given by 



ie-(o)e 

Jm,n\}J) — 



TT m! n! 



(55) 



Fi is the Gauss hypergeometric function [4a |. 70 is a parameter which depends upon the system-reservoir coupling 
strength and Nth is as defined above. Also 

-^ni — n (0 



n sinh(7o At/2) + A cosh(7o At/2) ' 



(56) 



5m — n(^) 



2Nth 



Q + Acoth(7oA</2)' 



2k. 



= Qm-n = 1 + 2Nth — (m - n), 
70 



and 



A = A„,-n = [n^ - mh{Nth + 1)] 



1/2 



(57) 
(58) 

(59) 
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p(e) 




e 



-3-2-10 1 2 3 



FIG. 6: Phase distribution Vifl) as a function of 6 for the dissipative anharmonic oscillator initially in a coherent state. The 
bold curve represents unitary evolution while the other curves represent temperature T = and evolution times f = 1,5 and 
10 for the dotted, large-dashed and dot-dashed curves, respectively. The parameters used are a) = 1, |C(0)P = 2, (fjo — Q (|C(0)|, 
cj>o are the initial state parameters), 70 = 0.01 and k = 0.05 (|52|) . 



H, R, X 




Ko I 



(a) 

FIG. 7: Number entropy H[m\ (large-dashed line), phase knowledge R[ff\ (small-dashed line) and entropy excess X[m,,9] (bold 
line) for an anharmonic oscillator system, initially in a coherent state, with respect to initial state parameter |^o| (= |?(0)|), 
subjected to a dissipative interaction where the parameters are as in the above figure. The figures depict an evolution time 
t = 2 and temperature T = 0. It can be shown that, in contrast to the corresponding harmonic oscillator case (Fig. [Ua)), 
increase in average number (|^o|) is not accompanied by a corresponding increase in phase knowledge. 



The corresponding complementary number distribution is obtained, using Eq. (j4]), as 

pirn) = Tmlfrn,Tn{t), (60) 

where /m,™ can be obtained from Eq. (|54l) . 

The phase distribution P{0) (|53|) is plotted in Figure |6] from which we see that with increase in time, phase gets 
randomized resulting in phase diffusion. Using this V{9) ([55]) in Eq. ^ to get the phase knowledge, p{m) (pO)) to 
get the number entropy and using these in Eq. (fTTj) we get the entropy excess which are plotted in Figure [T] The 
effect of the dissipative interaction is seen to manifest in the increased phase randomization and entropy excess with 
increase in the state parameter j^ol- The principle of entropy excess, Eq. pT|) . is, again, clearly seen to be satisfied. 



V. ATOMIC SYSTEM 



Here we study entropy excess (jlip for number-phase complementarity in (finite-level) atomic systems, briefly re- 
visiting results obtained from the perspective of an upper bound on the knowledge-sum of complementary variables 
in Ref. [l^ . An interesting generalization of the knowledge-sum of complementary variables could be made, in the 
context of quantum communication, using the information exclusion relations developed in j4l| . As pointed out ear- 
lier, the knowledge-sum approach cannot be applied to infinite dimensional systems, whereas the principle of entropy 
excess can be applied to finite as well as infinite dimensional systems, making it a more flexible tool for studying 
number-phase complementarity in a host of systems. 

As an application of the principle, it is appropriate to study the effect of noise. This we do for noise from both 
non-dissipativc as well as dissipative interactions of the atomic system S with its environment, which is modelled as 
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a bath of harmonic oscillators starting in a squeezed thermal state [T^, [10, [Hi • In Section fVBI we consider the effect 
of t he p hase damping channel, which is the information theoretic analogue of the non-dissipative open system effect 
[isl l5Cll |. while in Section IV CI we consider the effect of the squeezed generalized amplitude damping channel which is 
the information theoretic analogue of the dissipative open system effect [s^, [5l[ . 



A. The principle of entropy excess in atomic systems 



For a (noiseless) two- level (spin- 1/2) system, the plot of entropy H[m] for all atomic coherent states is given by 
the large-dashed curve in Figure[8l^a). The equatorial states on the Bloch sphere, corresponding to a' = tt/2, are the 
maximum knowledge state (MXK) states of (j), and are precisely equivalent to the minimum knowledge state (MNK) 
states of m (characterized by 7?[m] = 1), as can be seen from comparing the large-dashed and small-dashed curves 
in the Figure. Thus number and phase share with MUBs the reciprocal property that maximum knowledge of one 
of them is simultaneous with minimal knowledge of the other, but differs from MUBs in that the maximum possible 
knowledge of (f> is less than \og{d) = 1 bit, essentially on account of its POVM nature. 

Two variables form a quasi-MUB if any MXK state of either variable is an MNK state of the other, where the 
knowledge of the MXK state may be less than logc? bits. Thus, Jz and are quasi-MUB's (but not MUB's). 

From the dot-dashed curve in the Figures ([5]) , we numerically find an expression of the uncertainty principle to be 

X[m, <j)] EE H[m] - R[(t)] > (61) 

for all pure states in C^, in conformity with Eq. ^ and hence also in agreement with the principle of entropy excess 
([TT|). As (/) is a POVM but m represents a regular Hermitian observable, in general X[m, 0] 7^ X[(j>, m]. The inequality 
is saturated only for the Wigner-Dicke states (as seen from the dot-dashed curve in the Figure), when H[m] and R[<j>] 
identically vanish. 

As an expression of the uncertainty principle, the relation ((6T|) still leaves some room for improvement. First, it is 
not a tight bound. In particular, for equatorial states it permits R[(j)\ to be as high as 1, whereas as seen from the 
small-dashed curve in Figure [51 the maximum value of i?[(/>], which is w 0.245. 

Following Ref. [l^, one way to address this problem is to modify (pT|) to the inequality 

X^[m, 0] = H[m] - > (62) 

for all pure states in C^, where parameter ^ (> 0) is chosen to be the largest value such that inequality is satified 
over all state space. Through a numerical search, we found that /i « 4.085 for dimension d = 2 and ^ ~ 1.973 for 
d = 4. From the concavity of H[m] and the convexity of R[(j)], it follows that Eq. holds for any mixed state. 

The small-dashed and dotted curves are, respectively, and Comparing their corresponding curves in the 

Figure, we note the tighter bound imposed by X^[to, 0] than X[m, (j)]. 

B. Application to the phase damping channel 

The 'number' and phase distributions for a qubit, Hs = ^'^z: starting from an atomic coherent state \a' , 13'), and 
subjected to a phase damping channel due to its interaction with a squeezed thermal bath, are [isl [28| 



p{m) = (^^.^J^)(sin(a72))2(^+'")(cos(a72))2(^- 



Vicj)) = 7^ fl + T sin(a') cos(/3' + ujt - •^e-^^'^^''^^*)] . (63) 



2tt 



4 



For completeness, the function ^{t) appearing in the above equation is given in Appendix A. We note the symmetry 
preserved in Figures (|S]) (a) and (b), about a' = tt/2, the equatorial states. In the case of R[(j)], this is because of 
symmetry of sin(a'), as in Eq. (|63p for 7^(0) about 7r/2, whereas in the case of H[m\, the symmetry comes about 
because the cos(-) and sin(-) functions, in Eq. ()63|) for p{m), appear only as an even power. For QND interaction 
p{m) is time- invariant, whereas V{(f>) evolves in a way that does not affect this symmetry. 

Figure |Hl[b) depicts the effect of phase damping noise on the number entropy H[m] (obtained by using the number 
distribution p(m) (|63)) ). phase knowledge (obtained by using the phase distribution 7^(0) (|63|) ) , /^i?[(/>], X[m,(f>] 
(by using Eq. (|lip ) and X^[m, tj)]. Comparing it with the noiseless case, as in Figure[5Ja), we find that H[ni\ remains 
invariant because p{m) is not affected when a system undergoes a QND interaction, but there is an increase in X^ 
because of phase randomization with time. 
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FIG. 8: Entropy excess of a two level system subjected to QND interaction starting in an atomic coherent state |a',/3'), as 
a function of a', with /3' — 0.0. The large-dashed (resp., small-dashed) line represents H[m] (resp., i?[(?!>]). The dotted-curve 
represents i-iR[(j}] (where /i — 4.085). The solid (resp., dot-dashed) curve represents the entropy excess (resp. X); (a) depicts 
the noiseless case . There is no /3-dependence; (b) depicts the case of QND interaction. The parameters used are cj — 1.0, 
ojc ~ 100, 70 (|A1[) = 0.025, bath squeezing parameters (|A5|) r = 0.5 and a = 0. The plots in the figure (b) are for a temperature 
T = 10 and an evolution time t = 1. We note the symmetry in the figure about a' = tt/2. In (a), the points where = 0, 
namely the polar and the equatorial states, represent the coherent state. If complementarity is expressed in terms of knowledge 
sum [ij], these states correspond to maximum knowledge states. 



C. Application to the squeezed generalized amplitude damping channel 



The 'number' and phase distributions for a qubit starting from an atomic coherent state |a',/3'), and subjected to 
a squeezed generalized amplitude damping channel [5l| due to its interaction with a squeezed thermal bath, are [1^ , 



p(m = l/2,t) = i 



l_2^lUfl + ^)e-T"* 



sin2(a72) + I^(l-e-^'* 



(64) 



and 



1 

2^ 



H sin(a')<^ a cosh(Q;t) cos(0 — + uj smh( at) sin(0 — /?') 

4a L 

e 2 



(65) 



A derivation of Eqs. and (|55|) can be found in Ref. [231. For completeness, the parameters appearing in these 
equations are given in Appendix [B] 

Figures ini^a) and (b) depict the effect of squeezed generalized amplitude damping noise on the functions depicted in 
the noiseless case of Figure [SJ^a), without and with bath squeezing, respectively. Comparing them with the noiseless 
case, wc find as expected that noise impairs both number and phase knowledge. If the dependence on /?' is taken into 
consideration (cf. Ref. [1^), it can be shown that squeezing has the beneficial effect of relatively improving phase 
knowledge for certain regimes of the parameter space, and impairing them in others. This property can be shown to 



improve the classical channel capacity [51|. Further, bath squeezing can be shown to render i?[</>] dependent on /3'. 
On the other hand, it follows from Eq. (|64p that R[m] is independent of /3', so that X[m,(/)] is dependent on /3'. This 
stands in contrast to that of the phase damping channel, where inspite of squeezing, X[m,<p\ remains independent of 
/?' and, furthermore, squeezing impairs knowledge of (f> in all regimes of the parameter space. 



VI. DISCUSSIONS AND CONCLUSIONS 



In this paper, we have recast the number-phase complementarity for finite dimensional atomic as well as infinite 
dimensional oscillator, discrete (Hermitian) as well as continuous (positive operator) valued, systems as a lower bound 
on an entropic measure called the entropy excess. For maximally complementary systems, the bound is 0, independent 
of the system dimension. This is in contrast to the conventional entropy sum principle, which has a lower bound of 
log d. To tighten the constraint imposed by the bound on i? [</)], we replace this quantity by where /i is a positive 

number with values (approximately) 4, 2 and 1 for two-, four- and infinite-dimensional systems. Thus dimensional 
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a' 



(a) (b) 

FIG. 9: Entropy excess of a two level system subjected to a dissipative interaction starting in an atomic coherent state |a',/3'), 
as a function of a', with /3' = 0.0. The large-dashed (resp., small-dashed) line represents H[m] (resp., R[(f>\). The dotted-curve 
represents iiR[(:f>\. The solid (resp., dot-dashed) curve represents the entropy excess (resp. X). Here ui = 1.0, Uc — 100, <& 
HB5|) = 7r/8, 7o ~ 0.025, with the evolution time t — 1 and temperature T = 10: (a) bath squeezing parameter r (|B5|) =0: (b) 
r — 1. Comparison of (b) with (a) shows that squeezing impairs both number and phase knowledge, leading to an increase in 
the entropy excess (and X). 

dependence of the inequality enters indirectly through the form ^ = ^J.{d). Encouraged by the above numerical- 
analytical pattern, we conjecture that as the system dimension increases from two to infinity, jj falls monotonically 
from about 4 to 1. 

In this work, we have made precise the sense in which the variables p{m) and P{4>) may be thought of as or differ 
from conventional complementary variables (soj . There arc two main differences as follows. First: whereas states of 
maximum number knowledge (the eigenstates of the number operator) have the maximum knowledge of logd bits, 
the maximum phase knowledge states have less than logd bits, phase being a POVM. This was the motivation for 
introducing the weight quantity fi. Second: even more remarkably, states of maximum phase knowledge do not 
correspond to equal amplitude superpositions of number states. In other words, the unbiasedness is not mutual, but 
one-way, a situation we characterize as one-way unbiased bases (l^ . 

In the second aspect of our work, the above analysis is applied to physically relevant initial conditions of the system 
for unitary as well as non- unitary evolution, due to the interaction of the system with its environment. The system- 
reservoir interactions are chosen such that both dephasing (decoherence without dissipation) as well as dissipative 
(decoherence with dissipation) effects on the system evolution are studied. 

Some interesting features seen were, for e.g., a harmonic oscillator starting out from an initial superposition of 
coherent (cat) states. The entropy excess principle was seen to provide an interesting complcmentaristic characteriza- 
tion of the even and the odd cat states, in that the excess is almost zero for the even state, indicating that ignorance 
of number approximately equals phase knowledge while in the odd state, the entropy excess is finite indicating that 
there the ignorance of number significantly outweighs knowledge of phase. 

Our entropy-based formalism can modify current approaches to number-phase complementarity: e.g., one can study 
complementarity in conjunction with such phenomena as nonlinearity induced coherences and atomic squeezing in an 
effectively finite-level atomic system. In the conventional approach, complementarity can be graphically demonstrated 
by the constrasting behavior of the number and phase distributions (eg.. Figs 1 and 2 of Ref. [3(j|). 

As a concrete application to a finite dimensional system in an experimental scenario, we consider the energy 
manifold of the four levels of (for instance) ^^Rb atom. This is first mapped to a pseudo-spin system of spin 3/2 
while the effect of selection rules of atomic transitions in ^^Rb is preserved [H^]. Complementarity can then be 
studied using (entropic) knowledge of the number and phase variables as a function of laser detuning and vis-a-vis 
atomic phenomena like coherent population trapping (CPT) or elcctromagnetically induced transparency (KIT). For 
example, simulations indicate an increase in phase knowledge accompanying the formation of the CPT state. Noting 
that {6,(f>\J^\6,(j>) — j sin 0e"^ , where J_ = Jx — Jy, one can detect in a practical, intcrfcrometric set-up by applying 
J_ to one of the two interferometric arms implemented in an atom-laser system. With appropriate adjustments, (jj 
will then manifest as a phase shift in the interference pattern. 
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Appendix A: Some expressions pertaining to the phase damping channel 

For the case of an Ohmic bath with spectral density 

I{uj) = ^we-"/'^% 



(Al) 



where 70 and ujc are two bath parameters characterizing the quantum noise, it can shown that using Eq. (jAl|) one 
can obtain 03 



and 



-fit) = cosh(2r) ln(l + Lolt^) - — sinh(2r) In 
27r An 



V{t) = -— tan-\cj,t), 

(l + 4a;g(t-a)^) 
(l + c.2(t-2a)2f 



(A2) 



- ^ sinh(2r) ln(l + Aa^col), (A3) 

4:71 



in the T = hmit, where the resuhing integrals are defined only for t > 2a. In the high T limit, 7(t) can be shown 
to be in 



7(t) ~ — cosh(2r) 



2a;citan ^ (ujct) + \n ( ^-^) 



— smh(2rj 

2TrhuJr 



4:UJc{t ~ a)ta.n ^ {2ujc{t ~ aj) 



4:UJc{t ~ 2a) tan {ujc{t — 2a)) + 4aWc tan (2aWc) + In 



[1 +a;2(t- 2a)2" 
[l + 4a;2(t-a)2 



In 



1 



1 + 4a2w? 



, (A4) 



where, again, the resulting integrals are defined for t > 2a. Here we have for simplicity taken the squeezed bath 
parameters as 



cosh(2r(aj)) = cosh(2r), sinh (2r(a;)) = sinh(2r), 
^(lu) = aio, 



(A5) 



where a is a constant depending upon the squeezed bath. The results pertaining to a thermal bath can be obtained 
from the above equations by setting the squeezing parameters r and $ to zero. 



Appendix B: Some expressions pertaining to the squeezed generalized amplitude damping channel 

Here the reduced dynamics of the two level atomic system interacting with a squeezed thermal bath under a 
weak Born-Markov and rotating wave approximation is studied. This implies that here the system interacts with its 
environment via a non-QND interaction, i.e., [Hg, Hgj^] ^ such that along with a loss in phase information, energy 
dissipation also takes place. 

The parameter a (Eq. ([65]) ) is given by 



Further 



and 



where 



7'' = 7o(2iV+l), 

7- = 7o^^, 
N = Afth(cosh^(r) + sinh'^(r)) + sinh^(r), 



(Bl) 

(B2) 
(B3) 
(B4) 



M = -- sinh(2r)e**(27Vth + 1) = xe'*, 



(B5) 
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and 

^th = si ■ (B6) 

Here A^th is the Planck distribution giving the number of thermal photons at the frequency to and r, $ are squeezing 
parameters. The analogous case of a thermal bath without squeezing can be obtained from the above expressions 
by setting these squeezing parameters to zero. 70 is a constant typically denoting the system-environment coupling 
strength. 
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